Temperature Profiles in Hamiltonian Heat Conduction 
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We study heat transport in the context of Hamiltonian and related stochastic models with nearest-neighbor 
coupling, and derive a universal law for the temperature profiles of a large class of such models. This law con- 
tains a parameter a, and is linear only when a = 1. The value of a depends on energy-exchange mechanisms, 
including the range of motion of tracer particles and their times of flight. 
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In nonequilibrium physics, the Fourier law is an example of 
a simple phenomenological principle whose molecular origin 
is very hard to explain. Idealizing homogeneous thin rods and 
wires with uniform cross-section as 1-D objects, this law says 
that heat flux is proportional to temperature gradient times 
heat conductivity. Ever since Fourier's pioneering work (H, 
physicists have tried to derive this law from first principles. 
The current state of the art is summarized in the excellent re- 
views IHIilH which all point to the need for a deeper the- 
oretical understanding beyond the many existing models and 
simulations. In this Letter, we report on new results for certain 
types of Hamiltonian systems and their stochastic realizations. 

To study the Fourier law in a Hamiltonian context, the most 
common setting is that of a chain of identical units comprised 
of disks, plates, penduli, and the like, coupled with short range 
forces between them. At its two ends, the chain is coupled to 
mechanisms simulating heat baths maintained at two different 
temperatures. Because the components of the chain are identi- 
cal, one may expect heat conductivity to be constant along the 
chain, so that by Fourier's law, the temperature profile is lin- 
ear. This seems to be the predominant thinking behind much 
of the recent work on Fourier's Law, although nonlinear pro- 
files are known to occur in other contexts |5]. 

In this Letter we point out some simple and natural mecha- 
nisms that lead to various profiles - both linear and nonlinear 
- in concrete Hamiltonian models, and derive a universal law 
for the profiles of these and other systems. 

Summary of results: 

(A) We consider Hamiltonian models consisting of a chain 
of energy storing devices (ESD) that are fixed in place and 
coupled to each other. For such a setup, we show that the 
temperature profile can be linear or nonlinear depending on 
the nature of the coupling. More specifically, we assume that 
energy exchange in the system is mediated by tracers, which 
move from ESD to ESD redistributing energy according to 
the rules introduced in [6, 7]. We find that the profile is lin- 
ear if energy transfer is carried out by a single tracer that 
moves freely along the chain, whereas heat conductivity is 
temperature-dependent ( and hence the profile is nonlinear) if 
the tracers are confined to specific regions. 

(B) Our second result is a universal law that holds for very 
general coupled chains of Hamiltonian or stochastic systems 



with nearest-neighbor interaction (including those considered 
in (A)). We show that as the number of constituent cells goes 
to infinity, the stationary temperature profile is given by 
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(1) 



where 7l and Tr are the temperatures imposed at the left and 
right ends, x is the coordinate along the system (normalized to 
x £ [0, 1]), and a is a constant. The form of this law does not 
depend on details of the system (precise conditions are given 
later). The value of a in (1), however, depends on the nature 
of the coupling. 

In the case of locked-in tracers, i.e., tracers confined to the 
regions between neighboring ESD, the value of a depends on 
their time of flight. In particular, a = | if the energy is purely 
kinetic. In the case of a single tracer moving freely along the 
chain, a = 1 independent of its time of flight. We show the 
results above for stochastic realizations of the models in (A), 
and explain why one should expect these results to carry over 
to their Hamiltonian counterparts. 

Our primary concern in this Letter is the temperature profile 
(TP). The existence of local thermal equilibrium (LTE), which 
is important for the definition of temperature, will be proved 
in 1 8 ] in the infinite volume limit for one of the models treated 
here. (For results on LTE for other models, see i9l flolHlll .') 

Before embarking on specifics, we note again that in our 
models, we distinguish between communicating agents (CA), 
which in our case are tracer particles, and ESD, which in our 
case are turning disks with fixed centers (we emphasize that 
the ESD are not infinite reservoirs). In a real-world conductor, 
the difference between CA and ESD is often blurred, but for a 
good theoretical understanding it is useful to keep them sepa- 
rate. These concepts are distilled from the following beautiful 
model: 

The MLL Model 00]: We describe this model in some 
detail, as it contains the basic ingredients of the models in (A). 
The MLL Model is purely Hamiltonian, and very careful sim- 
ulations show that the Fourier law holds. The system consists 
of an arrangement of N disks of radius 1 placed as in FigQ 
and a little point particle of mass 1 (the tracer) which wan- 
ders around the playground £2 (the physical space occupied 
by the system minus the disks), bouncing off the disks 11711 . 
While FigQsuggests a Lorentz gas fl^l . there is a crucial dif- 
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FIG. 1: A typical arrangement of disks in the MLL Model. The sim- 
ulations in 01 were done with 2 rows and periodic boundary condi- 
tions in the vertical direction, and zig-zag reflecting walls of temper- 
ature 7l, resp. 7r, at the two ends. The tracer particle is not shown. 

ference here: Each disk is "nailed down" in its center, around 
which it turns freely. The state of the system is described by 
x = (a>i , . . . , CON,q, v) where to,- is the angular velocity of disk 
i, q^Q. is the position of the tracer, and v is its velocity. When 
the tracer collides with a disk, the rule of interaction is that of 
"sticky reflection": Suppose the angular velocity of the disk 
being hit is CO, and v n and v t are the normal, resp. tangential, 
component of v relative to the impact point. Then the values of 
v and CO after the collision are given by the energy and angular 
momentum conserving law 

V n = -V n , V t = Vt-y— (Vt-O)) , 

2 

to' = 03 + — — (vt-ffl) . (2) 

1 +£ 

Here, e G R + is proportional to the moment of inertia of the 
disk divided by the mass of the tracer and the square of the 
radius of the disk. |7] treats mostly the case £ = 1, where v[ = 
CO and co' = v t , i.e., the two quantities are simply exchanged. 
Of particular interest to us is the case e <C 1, which from the 
tracer's point of view resembles the classical Lorentz gas. 

Remark. The MLL Model, as well as the Hamiltonian mod- 
els we will describe later, have the following important prop- 
erty: Since there is only a hard-core potential, the time evolu- 
tion of the system is rescaled by VX when the energy of the 
particle and the disks are rescaled by A. In this respect, the 
model in |7] is very different from models such as the ding- 
a-ling and ding-dong models (HQ HI. Most importantly, 
the energies of the tracer and the disks alone determine the 
time-of-flight of the tracer: It does not depend on the history 
(as it would in many models considered so far 0QQ). 

We still need to say what happens when the tracer hits one 
of the ends. In |7|], many variants are considered, but for our 
purpose, the following process is assumed: When the tracer 
hits one of the ends, it exits the system, and a new tracer is 
injected into Q. to take its place. The new tracer enters Q. at 
the point of exit of the old one. Its direction is arbitrary, and its 
speed is given by the Maxwell distribution for the temperature 
of the end in question. 

Introducing the models studied in this Letter: We now 
introduce two classes of models that have the same basic setup 
of tracers and turning disks (and the same rules of interaction) 



as in the MLL Model. However, the configurations of disks 
and tracers in these two models are chosen to give rise to two 
conceptually very different modes of transport. 

Model I (Wandering tracer): A single tracer wanders 
along a chain of boxes separated by walls with a tiny hole that 
allows the tracer to pass between adjacent boxes. Deep inside 
each box is a turning disk surrounded by many fixed disks. 
The turning disk serves as ESD, while the fixed disks are bona 
fide Lorentz scatterers, which serve to randomize the angles of 
incidence in collisions between the tracer and the turning disk, 
leading to the exchange of a random portion (i.e., the tangen- 
tial component) of the energy of the tracer. The smallness of 
the holes in the separating walls keeps the tracer in each box 
for a long period. This together with the chaotic action of the 
Lorentz scatterers ensures that the tracer is equally likely to 
exit the box from either side. 

Stochastic realization: We model this system on a 1-D lat- 
tice with N sites. There is a random variable <!;,■ represent- 
ing the energy at site i, i= l,...,N, with values in [0,°°). 
The heat baths at the ends are modeled by stochastic vari- 
ables 4l an d <Hr which take values in [0,°°) with a distribu- 
tion 7Lexp(— %l/T]_), resp. TRexp(— ^r/Tr) (the Boltzmann 
constant Icq being set to 1). We identify with the variable 
£,0, and <^r with <^v+i- There are two more random variables, 
tj, to be thought of as the energy of the tracer, and i, which 
gives the location of the tracer at any given time. We assume 
that when the tracer is at site i, it interacts with At a site i, 
i $ {0,N + 1}, the action is as follows: There is a clock which 
rings with rate /, for example / = tj ~ 1 ' 2 or (77 + '/ 2 , rep- 
resenting the time it takes for the tracer to make its way around 
the ith box. When the clock rings, the following mixing of 
energies takes place 1 18]: Choose a random variable p with 
uniform distribution in [0,1]. Then, 

»?' =/>(&+ 17), U = (i-pm+n)- (3) 

If ; = 0, then tj is replaced by a value chosen from the ex- 
ponential distribution for the temperature T^. The rule at the 
right end (i = N+ 1) is similar. After these operations, the 
tracer jumps with probability j to i — 1 or i+ 1 - except when 
it is at the ends, in which case it stands still with probability 
5 or moves into position 1 (resp. N) from the boundary. We 
prove in [8] that for all such models, the TP is linear and the 
distributions of the satisfy LTE. 

Model II (Locked-in tracers): In this model there is a 
channel one-disk wide, with reflecting boundaries, and with 
turning disks located at fixed distances apart. These disks turn 
freely, but they block the channel completely, separating it 
into individual cells. Inside each cell is a single tracer, which 
moves back and forth, transferring energy between the turn- 
ing disks that border the cell (l^l . Here one can assume the 
tracer hits the two turning disks alternately, or, to further ran- 
domize the situation, one can add a number of fixed disks in 
each cell as illustrated in Fig.|2] After hitting one turning disk, 
the tracer then "gets lost" in this array of Lorentz scatterers, 
to emerge at some random moment to hit the turning disk at 
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FIG. 2: A sketch of Model II, when it is made more chaotic. Be- 
tween the rotating disks, there are disks serving as Lorentz scatterers. 
The tracers are not shown. The horizontal walls are reflecting. 



either side with equal probability. In both cases, the time-of- 
flight of the tracer between hitting turning disks depends only 
on the speed of the tracer (and not on the state of the disks). 

Stochastic realization: We have N sites on a 1 -D lattice as in 
the previous stochastic model, with representing the energy 
of the turning disk at site i. In this model, however, there is one 
independent variable 77, for each pair (/,/ +1), i = 0, . . . ,N, 
representing the energies of the tracers. Each site is equipped 
with an (independent) clock which rings at an exponential rate 

-1/2 

proportional to rj ; . When this clock rings, an exchange of 
energy involving 77, takes place. As in the Hamiltonian model, 
one may assume 77, exchanges energy alternately with ^, and 
or, rjj chooses with probability j its left or right partner 
(i.e., or ^,+i), and performs the usual mixing: For example, 
if £,i has been chosen, then 



tj,), # = (1-/0(6+71,) 



(4) 



When the clock at site rings, 770 is replaced by a value chosen 
from the exponential distribution of temperature 7l as before. 
The rule at the right end (i — N + 1) is similar. Numerical 
simulations show clearly profiles deviating strongly from lin- 
earity. They are in perfect agreement with the value of (X = | 
predicted by theory (see Fig|3]and the sketch of argument be- 
low). 

The qualitative shape of the nonlinear profile can be under- 
stood easily by considering 3 successive sites, say 
and Since <!;,'_ 1 < 1<+i rattles faster than 77,- (the 

rate being given by r\~ x l 2 ). Thus ^ equilibrates more often, 
and hence better, with than with which explains the 
concavity of the TP. 

We next present some details of the theoretical arguments 
for the results above [8]: 

Reduction from Hamiltonian to stochastic models: As 
dynamical systems, the Hamiltonian models above are very 
chaotic. With e in the MLL Model taken to be <C 1, the dy- 
namics are close to those in billiards of Sinai type and Lorentz 
gases. This chaotic behavior is used to induce a strong "mem- 
ory loss" for the tracer particles. Fast correlation decay fl^l 
ensures that the tracer "forgets" from which side it enters a 
site. This justifies our assumption in the stochastic model that 
the fraction of energy exchanged is independent from one step 
to the next. 

Range of validity of (1). We consider a coupled chain of 
length N ma. steady state, and let Ej be the mean temperature 
at site i. We assume (1) translation invariance of the model; 
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FIG. 3: The mean values of t, and rj, as a function of x = i/N, with 
Th = 1, ?k = 10, N = 100, averaged over 2 • 10 9 exchanges of energy 
according to Model II. Superposed is the theoretical curve of Eq.{TJ. 
Note that the temperature profile is not linear and that its curvature 
is more pronounced at the cold end. Here, a = | . We have obtained 
the same profile for many variants of Model II. 



(11) Ej is determined by Ej + \ and (111) energy-scale in- 
variance, meaning if £,±1 are multiplied by X, then Ej is also 
multiplied by A; (tv) as N — > °°, the TP tends to a smooth 
function. It is not hard to show that (i)-(iv) imply that the limit 
stationary density satisfies a second order differential equation 
(namely[6}, the solution of which with boundary conditions Ti 
and Is is (1). 

We demonstrate how to derive this law - and compute the 
constant a at the same time - for Model II: 

Computing the temperature profile in Model II. Con- 
sider 3 successive sites. For illustration we assume (i) each 
tracer visits alternately the left and right disks and (ii) the 
mixing of energies at each collision is exactly half-and-half. 
The stationarity condition means that the speed of the tracer is 
equilibrated as well. Let the mean energy of the left tracer be 
77 _ _> as it heads toward site i and 77 _ ^_ as it goes away from 

it. By the rule of mixing, we have 77 » = (77_.^_ +2?j_i), 

TJ__ = (t?_^ +E-), leading to T7 J = (2eU + Et)/3, 
tj_«- = (£,-_! +2£,-)/3. 

We assume the speed of the tracer is E" 1 when E is its en- 
ergy. The value of 7 for Model II as described above (and 
many other models without potential) is 7 = 5 (since the en- 
ergy is purely kinetic) while for potential interactions, the time 
is given by an integral of the form f dq (E — V (q r )) -1 / 2 , which 
for large E and V(q) ~ \q\'" behaves like (?(E~ l l 2+l l m ), so 
that 7 = i — — . Fixing 7, the average time for a round-trip 



of the tracer between sites i — 1 and i is T_ = 77 



-7 



and the rate at which site i gets information from the left is 
the inverse of this quantity. An entirely analogous reasoning 
applies on the "+" side. From the stationarity condition, we 
get 

p TI 1 (£, + 77_,^)/2 + Tf 1 {Ej + 77 + ,^)/2 

t{ = j = . p) 

T_ + 1+ 

Perform a perturbative analysis at the point x = i/N where 
is very large. Then, to second order in e = l/N, Ei±\ = 
T(x)±eT'{x) + \e 2 T"(x) , and (0 leads to 



t = T(x) + e : 



T"{x)T(x) + r(T'(x)Y 
4T{x) 



+ 0(e 4 ) . 
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Since £, must be equal to T(x), we find 

T"(x)T(x) = -y{T\x)f , (6) 

the solution of which with boundary conditions T(0) = 7T and 
T(l) = T R is Eq.Q with a = 1 + y. Thus for Model II, y = \ 
and a = I. One also checks that the energy flux is given 
by T' \x)^/t(x) (which is constant along the profile, but not 
proportional to the temperature difference). 

Generalization. Note that when 7=0, i.e., when the rate 
at which information is exchanged is independent of energy, 
then a = 1, which indeed gives a linear TP. Note also that our 
derivation is quite general: if E r is replaced by 1 /F(E), the 
profile is given by T"(x)F(T(x)) = (T'(x)) 2 F'(T(x)). 

Remark. Many authors have done careful simulations of 
models that are close to Model II, and have observed linear 
TPs. It should be noted that the profiles predicted by Q are 
very close to linear when Tr/Tl is not far from 1. Our theory 
predicts that deviations from linearity become more promi- 
nent with the increase of Tr/T^. 

Linear profile and LTE for Model L The reasoning above 
is not valid for Model I, for here it is a single tracer that is 
responsible for all transmission of information to all sites. 

In the discussion below, we refer to the stochastic realiza- 
tion of Model I as a variable-time model, in the sense that the 
tracer lingers in each site for variable time periods depend- 
ing on the function /. We introduce an associated fixed-time 
model, obtained by neglecting the waiting time above, i.e., all 
rules are as in the variable-time model except that the tracer 
jumps at fixed time intervals. 

Linear profile. Consider first the fixed-time model. Here 
the tracer performs a standard random walk on {1,2, .. . ,N} 
(except at the ends). It is easy to see that it spends an equal 
amount of time at each of the N sites. Indeed for a fixed site 
i, the tracer comes from the left and the right exactly the same 
number of times over any period, for if the tracer heads left 
from site i, then it can only return from the left; and similarly 
for the right. Reasoning as above, we see that in the stationary 
measure, is affected equally by E, r_ i and £ 1+ 1 , hence the TP 
is linear, i.e., a = 1. Similar reasoning applies to the variable- 
time model. Note that the mixing rule, and hence the outcome 
of the mixing process, does not depend on the waiting time. 

LTE. Methods similar to those in 11(1 ll ill with a different 
dual XXX give LTE for the fixed-time model. While there is 
a correspondence between the sample paths for the variable- 
time model and its fixed-time counterpart, the varying du- 
rations the tracer spends at each site affects nontrivially the 
stationary measure. Proving canonical distributions for the 
model of interest, i.e., the variable-time model, requires fur- 
ther control of correlations 1 8]. Our proof applies to any func- 
tion / that is integrable and local, meaning it depends only on 
the variables and 77 when the tracer is at site i (which is the 
case here). 

Remark: Even though this LTE result is for the stochastic 
realization of Model I, we have incorporated into our setting 



the notions of trajectories and velocities of particles, two of 
the key ingredients in Hamiltonian dynamics. 

Conclusion. ( 1 ) We have pinpointed some simple and very 
natural mechanisms responsible for both linear and nonlinear 
profiles in homogeneous conductors modelled by Hamiltonian 
systems. Our results show in particular that when their ranges 
of motion are restricted, interacting particles (or springs) ac- 
quire energy-dependent speeds. (2) We have derived an exact 
formula - a universal law - for the energy profiles of very gen- 
eral Hamiltonian and stochastic chains with nearest-neighbor 
interactions. When nonlinear, this law predicts how deviation 
from linearity increases with the quotient 7r/7l. (3) Finally, 
the underlying causes for nonlinearity that we have identified 
clearly go beyond the models studied here. They suggest that 
the presence of some (weak) nonlinear effect may be a more 
common phenomenon than recognized when very disparate 
temperatures are imposed at the two ends of a 1-D system. 
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